Dust Static Spherically Symmetric 
Solution in f{R) Gravity 

Muhammad Sharif *and Hafiza Rizwana Kausar ^ 
Department of Mathematics, University of the Punjab, 
Quaid-e-Azam Campus, Lahore-54590, Pakistan. 

Abstract 

In this paper, we take dust matter and investigate static spherically 
symmetric solution of the field equations in metric f{R) gravity. The 
solution is found with constant Ricci scalar curvature and its energy 
distribution is evaluated by using Landau-Lifshitz energy-momentum 
complex. We also discuss the stability condition and constant scalar 
curvature condition for some specific popular choices of f{R) models 
in addition to their energy distribution. 
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1 Introduction 

Our latest data from different sources, such as Cosmic Microwave Background 
Radiations (CMBR) and supernova survey indicates that energy composition 
of the universe is the following: 4% ordinary matter, 20% dark matter and 
76% dark energy ^\ The dark energy has large negative pressure, while the 
pressure of the dark matter is negligible. The current accelerated expansion 
of the universe may be due to the presence of dark energy so-called effec- 
tive cosmological constant. There are various directions aimed to construct 
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the acceptable dark energy model. For example, quintessence or phantom 
models, dark fluid with complicated equation of state. String or M-thcory, 
higher dimensions, more complicated fleld theories, etc. Despite the number 
of attempts, there is still no satisfactory explanation about the origin of dark 
energy. 

The f{R) theory of gravity provides the very natural gravitational alter- 
native for dark energy. The cosmic acceleration can be directly explained by 
taking any negative power of the curvature (like term) In this way, this 
theory helps in modification of the model to achieve the consistency with the 
experimental tests of solar system. A unification of the early time infiation 
and late time acceleration is allowed in f{R) theory. It is also very useful 
in high energy physics for explaining the hierarchy problem and unification 
of GUTs with gravity ^\ 

This theory has produced a great number of papers in recent years, for ex- 
ample. Several features including solar system test ^\ Newtonian limit 
^\ gravitational stabihty ^'^^ and singularity problems are exhaustively 
discussed. Much work has been devoted to place constraints on f{R) models 
using the observation of CMB anisotropics and galaxy power spectrum 
Kobayashi and Maeda have studied relativistic stars in this theory grav- 
ity. It is shown explicitly that stars with strong gravitational fields develop 
curvature singularity and hence are prohibited. Erickcek et al. found 
the unique exterior solution for a stellar object by matching it with interior 
solution in the presence of matter sources. Kainulainen et al. ^^-^ studied the 
interior spacetime of stars in Palatini f{R) gravity. 

The vacuum solutions of the field equations in metric f{R) gravity has 
attracted many people. Since the spherically symmetry plays a fundamental 
role in understanding the nature of gravity, most of the solutions are discussed 
in this context. Multamaki and Vilja investigated static spherically sym- 
metric vacuum solutions of the field equations. It is shown that solution 
with constant scalar curvature corresponds to Schwarzschild de Sitter space- 
time for a specific choice of constants of integration. Carames and Bezerra 

discussed spherically symmetric vacuum solutions in higher dimensions. 
Capozziello et al. analyzed spherically symmetric solution using Noether 
symmetry. 

Azadi et al. have studied cyhndrically symmetric solutions in Weyl 
coordinates. They have shown that constant curvature solutions reduce to 
only one member of the Tian family in General Relativity (GR). Momeni 

has found that the exact constant scalar curvature solution in cylindri- 
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cal symmetry is applicable to the exterior of a string. In a recent work, 
Sharif and Shamir have studied exact solutions of Bianchi types / and V 
spacetimes in f{R) theory of gravity. Multamaki and Vilja investigated 
non-vacuum solutions by taking perfect fluid. It is found that for a given 
matter distribution and equation of state, one cannot determine the function 
form of f{R). HoUestein and Lobo explored exact solutions of the field 
equations coupled to non-linear electrodynamics. 

The energy localization in GR is a serious and long standing issue but 
without any definite answer. The well-known energy-momentum prescrip- 
tions are given by Landau-Lifshitz, Einstein, Bergmann, Papapetrou, Gold- 
berg and M0ller. Virbhadra et al. showed that five different energy- 
momentum complexes yield the same energy distribution for any Kerr-Schild 
class metric. Chang et al. proved that every energy-momentum complex is 
associated with a Hamiltonian boundary term Recently, this problem 
has been attempted in alternative theories of gravity. In this connection, 
reasonable amount of work has been pubhshed in teleparallel theory of 
gravity. Multamaki et al. are the pioneers to discuss energy problem in 
/(i?) theory of gravity. They presented the concept of energy-momentum 
complex (EMG) in this theory and generalized the Landau-Lifshitz energy- 
momentum complex. The prescription is used to evaluate energy-momentum 
for the Schwarzschild-de-Sitter spacetime. Recently, Sharif and Shamir 
found energy densities of plane symmetric and cosmic string spacetimes. 
They also discussed the stability conditions. 

The purpose of this paper is two fold: Firstly, we study non-vacuum 
static spherically symmetric solutions of the field equations using metric f{R) 
gravity in the presence of dust fiuid. Secondly, we use generalized Landau- 
Lifshitz energy-momentum complex to evaluate energy density for constant 
scalar curvature solution. The paper is organized as follows. In section 2, 
we present spherically symmetric field equations and some of the relevant 
quantities. Section 3 is devoted to study the non-trivial solution of the field 
equations. In section 4, we calculate the generalized Landau-Lifshitz energy- 
momentum complex for constant scalar curvature solution and also discuss 
some well-known /(/?) models in this context. In the last section 5, we 
summarize and discuss the results. 
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2 Field Equations in f{R) Gravity 



The action in f{R) gravity is given by 

'f{R) 



S = d Xy/^ 



2k 



+ L 



M 



(2.1) 



where f{R) is a function of the Ricci scalar and Lm is the matter Lagrangian 
depending upon the metric (^^j, and the matter fields. Variation of this action 
with respect to the metric tensor leads to the following fourth order partial 
differential equations 

F{R)R^, - ]^f{R)g^, - V^V,F{R) + g^,UF{R) = kT^,, (2.2) 

where F{R) = df{R)/dR, □ = V^V^ with representing the covariant 
derivative and k,{= 8tt) is the coupling constant in gravitational units. Taking 
trace of the above equation, we obtain 

F{R)R-2f{R) + 3nF{R) = 8itT. (2.3) 

Here R and T are related differentially and not algebraically as in GR {R = 
—kT). This indicates that the field equations of f{R) gravity will admit 
a larger variety of solutions than does GR. Further, T = does no longer 
implies i? = in this theory. 

The Ricci scalar curvature function f{R) can be expressed in terms of its 
derivative as follows 

/(«) = -»-T + F(R)R + 3nF(R)^ ^ 

This is used to study various aspects of f{R) gravity, particularly its stability, 
weak field limit etc. Substituting this value of f{R) in Eq. fl2.2p . we obtain 



FiR)R - DFjR) - 8nT _ F{R)R^^ - V ^F{R) - SttT^^ 

4 gt^t. 



(2.5) 



In the above equation, the expression on the left hand side is independent of 
the index /i, so the field equations can be written as 

^ _ F{R)R^^-V^V^F{R)-87rT^, ^ ^^.e) 
Notice that is not a 4- vector rather just a notation for the traced quantity. 
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2.1 Spherically symmetric spacetime 

We take the following static spherically symmetric spacetime 

ds^ = A{r)dt^ - B{r)dr'^ - r\de^ + sin^ ed(j)'^). (2.7) 

The components of the Ricci tensor are 

1 4'2 A'R' A' 

1 , A'B' B'A^ , , 

i?33 = sin2^i?22. (2.11) 

The corresponding Ricci scalar is 

-2, r'^A" A', PA\ B' rM\, 

R=^—\l-B + + — r r + , 2.12 

r^B^ 2A A^ AA' B^ AA ^ ' 

where prime denotes derivative with respect to the radial coordinate r. The 
dust energy-momentum tensor is given as 

Ty^y = pUf,Uu, (2.13) 

where = 5° is the four-velocity in co-moving coordinates and p is the 
density. Since Eg. (12. 5 p is independent of index p, so A^ — A^ = for all p, v 
and yields the following two independent equations 

F" F' A' B\ A' B' Stto 

l-O D /I D«,J I^O/ID A Al ID A A TD2 n A TD^ 



2BA Br '2AB AA^B AAB^ 2ABr 2B^r 
1 1 Svrp 

r"2 ~ R^^ ~ 'a' 



+i^-iz:^]--r = o. (2.15) 



Thus we get a system of two non-linear differential equations with four un- 
known functions, namely, F{r), p{r), A{r) and B{r). 
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3 Solution of the Field Equations 



This section is devoted to discuss solution of the field equations by assuming 
constant scalar curvature which are directly involved in explaining the ac- 
celerating universe. In order to take into account acceleration of the present 
universe, we need to take very small value of the constant to f{R) 

The conservation law of energy-momentum tensor, T^.^, = 0, for dust 
matter gives A = constant = Aq (say). Thus the system of field equations 
f l2.14p and (12.151) is reduced to three unknowns with the following two non- 
linear differential equations 

-1 „ B' , B' 8np 

Using the assumption of constant scalar curvature {R = Rq), i.e., F{Rq) = 
constant, the field equations become 

|^FW-^=0. (3,18) 

' +i-iWi^)-5j^! = 0. (3.19) 



'2B^r r2 Br^' ' Ao 

Now we have two differential equations with two unknowns, B{r) and p(r). 
Subtracting Eq. (l3.19p from Eq. (l3.18l) . we have an ordinary differential equa- 
tion in terms of B{r) as follows 

B'r + 2B-2B^ = (3.20) 

which has the following solution 

Bir) = (3.21) 
1 — cir^ 

where Ci is a constant. Inserting this value of B in any of the above equations, 
it follows that 

P = ^ = Po (3.22) 

which is purely a constant. We would like to mention here that if matter den- 
sity and scalar curvature are constant and also po depends on gravitational 
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constant and effective cosmological constant then the field equations be- 
come equivalent to the Einstein field equations. The spacetime for constant 
curvature solution takes the following form 



ds^ = Aodt 



^ - . ^ J r^ - r\de^ + sin^ 9d(j)^). 



1 — Cir 



(3.23) 



This solution corresponds to the well-known Tolman-Oppenheimer-Volkoff 
(TOV) spacetime when density is constant and pressure is neglected For 
all ordinary, non-relativistic stars where p << p (for example at the center of 
the sun), we can neglect pressure and can consider an idealized object with 
a constant density. 

The scalar curvature turns out to be Rq = 6ci and hence Eq.f l2.3p yields 



For the acceptability of any f{R) model, it is necessary to satisfy this equa- 
tion. In the next section, we check this condition for some well-known f{R) 
models and then calculate energy density using the generalized Landau- 
Lifshitz energy-momentum complex. 

In view of the above information, the universe could start from inflation 
driven by the effective cosmological constant at the early stage where cur- 
vature is very large. With the passage of time, the effective cosmological 
constant also becomes smaller corresponding to the smaller curvature. After 
that time, the density of matter or radiations become small and curvature 
goes to constant value Rq. Thus expansion could start and cosmological 
constant can be identified as f{Ro) in the present accelerating era. 



4 Landau-Lifshitz Energy-Momentum Com- 
plex 



Now we evaluate energy density for the constant scalar curvature solution. 
For this purpose, we use the generalized Landau-Lifshitz energy-momentum 




(3.24) 



Substituting the values of po and Rq, it follows that 



/(i?o) = 2ci/'(i?o). 



(3.25) 
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complex. We would like to mention here that this energy-momentum complex 
is valid only for those solutions which have constant scalar curvature. 

The generalized Landau-Lifshitz energy-momentum complex is given 
as follows 

T>^^ = f\RoKl + ^(/'(^o)i?o - /(i?o))^(^7'^V - g^\n. (4.26) 

Here rj^^ is the Landau-Lifshitz energy-momentum complex evaluated in the 
framework of GR and is given as 

<l = {-9){T^'' + t>rL). (4.27) 
where t^'^ can be obtained by the following formula 



+ 9'''9'\r^.,ns-K,^',s)]- (4.28) 
The 00-component of Eq. f l4.26p yields 

OK OX^ 

= f{R,)rll + -{f'{Ro)Ro-f{Ro)){^x' + Sy'') (4.29) 

with 

rll = {-y){T'' + tT,). (4.30) 
We can find t° l from Eq. (14.281) as follows 



~ 2k '^'^ ^05^-/37 ^a^^l3S){9 9 9 9 ) 



I ^Oa /37|'p0 p<5 _|_ pO p(5 _ pO ^5 _ pO -nS \ 

+ y y l-"- a.S-'- ^7 + I3j^ aS ^ -yS^ af3 ^ al3^ -yS) 

+ fl' l-L a^J- ^7 + J- /37J- a<5 J- 7<5J- J- a/3-'- 75J 

+ y"^g^\r%r% - r%r?^s)], «, 7, ^ = 0, 1, 2, 3. 
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which finally gives 



tfL = ^[(l-ciO{^^ + ^^cot^-2r^ + icot^} 
IdttAo 1 — Cir^ 1 — cir^ r 



+ -4{ ^""^ - -cotg-2(l-cir^) -2cot'g} 
1 — cir^ r 



1 r 2ci 4 6 , 2ci 
r"' sm y 1 — cir"^ ^ ]^ _ (..j^^^ 

+ 2cot2^-2cos^^}]. (4.31) 

Inserting this value in Eq.f l4.30p and then substituting the resulting equation 
in Eq.(i29]), it follows that 

00 .//D N^O^^sin^^ p 1 2ci 



1 — cir^ IGttAo 1 — Cir^ 

'^Cl ^ 9 1 1 r 2ci 1 

+ ^ ^cot^-2r2 + -cote} + -{- ^--cot^ 



± 2ci 4 6 

2(1 - cir^) - 2 cot^ + . {- ^ + + - cot ^ 



+ cot ^ + 2 cot^ 6 - 2 cos^ 6}}] + —L^{f'{R,)R, 

- /(i?o)). (4.32) 

This is the energy density satisfying the condition of constant scalar curva- 
ture. We can evaluate this quantity for different well-known f{R) models. 
Also, we check the validity and the stability condition for these models in 
the context of cosmology. 

4.1 Energy Density of the First Model 

First of all we evaluate energy density for the model f{R) = R + eR^, where e 
is any positive real number. The stability criteria for this model is restricted 
to e < which corresponds to / (R) > 0. For e = 0, GR is recovered 
in which black holes are stable classically but not quantum mechanically 
due to Hawking radiations. Since such features also found in f{R) gravity, 
hence the classical stability condition for the Schwarzschild black hole can be 
enunciated as / [R) > 0. It is interesting to mention here that this model 
satisfies the condition of constant scalar curvature for a specific value of the 
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constant Ci, i.e., Ci = ^. Further, the stabihty condition for this model 

,1 „ , = ^ > 0, is also satisfied. 

The 00-component of the generalized EMC takes the form 



, 2ci rci ^ ^ 9 1 ^-1 1 r 2ci 

X {- -^ + - ^cot^-2r' + -cot^} + -{— ' 

1 — Cir"^ 1 — Cir^ r 1 — 



_ lcote-2(l -cir^) - 2 cot^e} + — ^{-^^ + 4 



6 2ci 9ec^ 
+ -cot^ + - ^cot^ + 2cot^^-2cos^^}}] + ^ 



r 1 — cir^ 477^40 

4.2 Energy Density of the Second Model 

Here we use the model f{R) = R — ^ — hE? to evaluate energy density, where 
a and h are any real numbers. This model involves the negative power of the 
curvature which supports the cosmic acceleration. In this way, any negative 
power of the curvature can be taken into account to achieve the consistency 
with experimental tests of Newtonian gravity. However, the model involving 
such term might not satisfy the stability conditions which can be significantly 
improved by adding square terms of the scalar curvature. For R — Rq — 6ci, 
we have 

f(Ro) = 6ci - ^ - 36bcl f{Ro) = 1 + ^ - 126ci. (4.33) 

The constant curvature condition implies that both the parameters a and b 
are related by the following expression 

1 - - Sbci = 0. (4.34) 

Further, the stability condition /"(Ro) > yields a + b{Ro)^ > which is 
satisfied for c\>^. Thus the model is acceptable for such choice of ci. 
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Consequently, the energy density takes the form 



2ci rci .1^91 
+ cote -2r^ + - 



+ ^{^^--cot^-2(l-cir')-2cot2^} 
1 — civ^ r 



1 r 2ci 4 6 , 2ci 
sm 1 — Cir^ ^2 ^ ]^ _ ^.^^^^ 

+ 2cot2^-2cos2^}}] + — ^(-^ -366c2). (4.35) 

IdttAq 3ci 

4.3 Energy Density of the Third Model 

The model considered here is f{R) = R - (-1)''-^^ + {-l)"'-^bR"', where 
m and n are positive integers and a, b are real numbers, and is widely used 
in cosmology. A model of such type with m = 1 — | with a depending upon 
the mass of the galaxy is used to approximate galaxies by taking spherically 
symmetric solutions. It is straightforward that one cannot fit the data for all 
astronomical masses for a single choice of f{R) as a depends upon the mass 
of individual galaxy. For R — Rq, the model becomes 

^ (6ci)"-^- - i-iy-'a + (-l)-^K6r,)"'+" 

(6ci)"+^ + (-l)"-^na + (-l)--^&m(6ci)"-+" 
= ^^^^y;^^^ . (4.37) 

For this form of /(i?), the constant curvature condition can be written as 

2(6ci)"+^ - (-l)"-'a(3 + n) + (-l)"-^6(3 - m)(6ci)"^+'^ = 0. (4.38) 
In particular, when m = 3 or 6 = 0, we get 

(-l)"-i2(6ci)"+^ 



n + 3 



n ^ 3. (4.39) 



It is worthwhile to mention here that the stability condition is satisfied for 
this value of a. 
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The corresponding 00-component gives 

(6ci)''+i + (-l)"^ina + (-l)™^i6m(6ci)"+" 



(6ci) 



n+l 



2ci rci ^ 9 1 ^, 

X {- ^ + 1 ^cot^-2r2 + -cot^} 

1 — cir^ 1 — cir^ r 

12 1 

+ — { ^--cote-2(l-cir^)-2cot^e} 

r 1 — Cir r 

1 , 2ci 4 6 , 2ci 
+ o • 2/, {l ^ + ^ + -cot^+:j ^cot^ 

+ 2cot^g-2cos^g}}]+ , 

^ i-ir-Mn + 1) + - l)(6ci)"-+" 

(6ci)" ■ ^ ■ ^ 

4.4 Energy Density of the Fourth Model 

This is an interesting model due to its logarithmic dependence on curvature 
and it also satisfies the existence of relativistic stars. It is given as follows, 
f{R) — R — aln(^) + (— l)"*~^6it!"^, where m is a positive integer, /c is a 
positive real number and a is any real number. For constant curvature Rq, 
we have 

fir 

f{R) = 6ci - aln(^) + (-l)™-^6(6ci)"^ (4.41) 

and 

^ 6c,-a+{-ir-^bm{6c,r _ 
6ci 

The constant curvature condition is expressed as 

4ci - a(ln(^) - 1) + (-1)-^6(1 - j){6c,r = 0. (4.43) 



For m = 3 or 6 = 0, this reduces to 



(4.44) 
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The corresponding 00-component is of the form 

6ci - a + (-l)'"-^6m(6ci)'"-i 



6ci 

X { ^ + ^cot^-2r^ + -cote} 



12 1 
+ — { ^--cot^-2(l-cir^)-2cot^^} 

r 1 — cir r 



1 r 2ci 4 6, 2ci 
sm 9 1 — Cir^ j.^ j- i _ ^^j.^ 

+ 2cot^^-2cos^^}}] + ^ 



X {a(ln(^)-l) + (-ir-^6(m-l)(6cir}. (4.39) 
k 

5 Outlook 

This paper investigates solution of static spherically symmetric spacetime 
with non-trivial matter distribution. We have restricted our analysis to the 
dust case and obtain solution with assumption of constant scalar curvature. 
In addition, we have explored energy localization problem using the general- 
ized Landau-Lifshitz energy- momentum complex in f{R) gravity. 

The scalar curvature for this solution turns out to be non-zero constant. 
This leads to constant density of dust matter and corresponds to the well- 
known Tolman-Oppenheimer-Volkoff spacetime when density is constant and 
pressure is neglected Thus it is interesting to discuss such solutions for 
idealized objects where density is constant and pressure can be neglected. 

For such solutions with constant curvature, we can use the generalized 
Landau-Lifshitz energy-momentum complex to discuss energy-momentum 
distribution. We have evaluated energy density for this solution as well as 
for certain specific f{R) models. It is worthwhile to mention here that the 
well-known f{R) models satisfy the stability condition as well as constant 
scalar curvature condition for this solution. The cosmological importance of 
all these models is also discussed. The model with negative power of the 
scalar curvature directly supports the cosmic acceleration. If we choose the 
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model with positive powers of the scalar curvature (higher than 1), the mass 
of the scalar field can be adjusted to be very large and the stability can be 
improved. For solutions with zero scalar curvature, the generalized Landau- 
Lifshitz EMC coincides with Landau-Lifshitz EMC in GR. This work adds 
some knowledge about the longstanding and crucial problem of the localiza- 
tion of energy. It gives the energy density expressions for important f{R) 
models which may help at some stage to overcome the theoretical difficulties 
in the cosmological and astrophysical context. 

The applicability of solutions could be tested by comparing with con- 
straints of the solar system and cosmology. It would be interesting to inves- 
tigate solutions for non-static spacetimes with energy-momentum tensor of 
other types of fluids. 
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